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Abstract. The discrete Laplace operator on a triangulated polyhedral surface 
is related to geometric properties of the surface. This paper studies extremum 
problems for eigenvalues of the discrete Laplace operators. Among all triangles, 
an equilateral triangle has the maximal first positive eigenvalue. Among all 
cyclic quadrilateral, a square has the maximal first positive eigenvalue. Among 
all cyclic n-gons, a regular one has the minimal value of the sum of all nontrivial 
eigenvalues and the minimal value of the product of all nontrivial eigenvalues. 



1. Introduction 

1.1. A polyhedral surface S* is a surface obtained by gluing Euclidean triangles. It 
is associated with a triangulation T. We assume that T is simplicial. Suppose (E, T) 
is a polyhedral surface so that V, E, F are sets of all vertices, edges and triangles 
in T. We identify vertices of T with indices, edges of T with pairs of indices and 
triangles of T with triples of indices. This means V ~ {1,2, ...\V\}, E = {ij \ i,j S 
V} and F = {Aijk \ i,j,k G V}. A vector (/i,/2, ...,/n/|)* indexed by the set of 
vertices V defines a piecewise-linear function over {S, T) by linear extension. 
The Dirichlet energy of a function / on is 



Es{f) = \ Jj^fl'dA. 



IS 

When / is obtained by linear extension of (/i, /2, the Dirichlet energy 
of / turns out to be 

%.T)(/) = ^ E - A)' + cot«L(A - + cot4.(/, - f,f] 

ijkeF 

where the sum runs over all triangles of T and for a triangle ijk S F, a'jf., aj.,^, a^j 
are angles opposite to the edges jk,ki, ij respectively. 

Collecting the terms in the sum above according to edges, we obtain 

ijeE 

where the sum runs over all edges of T and 

i(cot a^j + cot a\j) if ij is shared by two triangles ijk and ijl 
^{cot afj) if ij is contained in one triangles ijk 

The Dirichlet energy of a piecewise linear function on a polyhedral surface was 
introduced and the formula ([T]) was derived by R. J. Duffin [D], G. Dziuk [Dzj 
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and U. Pinkall & K. Polthier [PPj in different context. For application of the 
Dirichlet energy and formula ([T]) in the characterization of Delaunay triangulations, 
see [R1 [G1 IBS[ ICXGL] . For interesting application of the Dirichlet energy and 
formula ^ in computer graphics, see the survey [BKPAL) . 

1.2. The discrete Laplace operator L can be introduced by rewriting the Dirichlet 
energy using notation of matrices 



By definition L is positive semi-definite. Its eigenvalues are denoted by 

= Ao < Ai < A2 < ... < A|y|_i. 

The discrete Laplace operator and its eigenvalues are related to the geometric 
properties of the polyhedral surface (5, T). For example, it is proved in jCXGLj that 
among all triangulations, the Delaunay triangulation has the minimal eigenvalues. 
In [GGLZ| . it is shown that a polyhedral metric on a surface is determined up to 
scaling by its discrete Laplace operator. 

1.3. In smooth case, the spectral geometry is to relate geometric properties of a 
Riemannian manifold to the spectra of the Laplace operator on the manifold. One 
of the interesting result is the following one due to G. Polya. For reference, for 
example, see [H], page 50. 

Theorem (Polya). The equilateral triangle has the least first eigenvalue among 
all triangles of given area. The square has the least first eigenvalue among all 
quadrilaterals of given area. 

It is conjectured that, for n > 5, the regular n-gon has the least first eigenvalue 
among all n-gons of given area. 

1.4. In this paper, similar results as Polya's theorem are obtained for the discrete 
Laplace operator. 

Theorem 1. Among all triangles, an equilateral triangle has the maximal Ai, the 
minimal A2 and the minimal Ai + A2. 

A cyclic polygon is a polygon whose vertices are on a common circle. By adding 
diagonals, a cyclic polygon is decomposed into a union of triangles. For each inner 
edge of any triangulation of a a cyclic polygon, the weight Wij is zero. Therefore 
the discrete Laplace operator is independent of the choice of a triangulation of a 
cyclic polygon. 

Theorem 2. Among all cyclic quadrilaterals, a square has the maximal Xi, the 
minimal Ai + A2 + A3, the minimal A1A2 + A2A3 + A3A1 and the minimal A1A2A3. 

Theorem 3. For n > 5, among all cyclic n-gons, a regular n-gon has the minimal 
^"^^^ Xi and the minimal 11"=!^ ^i- 



E^S,T){f)^^{fl,-,f\V\Wl,...,fv\) 



where each entry of the matrix L is given as 
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1.5. Plan of the paper. Theorem [l] Theorem [2] and Theorem |3] are proved m 
section 2, section 3 and section 4 respectively. 

2. Triangles 

2.1. In this section we prove Theorem [TJ Let 61,62,63 be the three angles of a 
triangle. Let a,; := cot 6i for i = 1, 2, 3. The condition 61 + 62 + 63 = tt implies that 

(2) 0102 + 0203 + 0301 = 1. 

The discrete Laplace operator is 

(Oi + O3 -Oi -O3 

-Oi Oi + 02 -02 
-03 -02 02 + 03 

The characteristic polynomial of Lg, is 

Psix) — det(L3 — XI3) — —x^ + 2(oi + 02 + 03)0;^ — 3(oi02 + 02O3 + 0301)0; 
= —x^ + 2(oi + 02 + 03)2;^ — 3a; 

by the equation 

The eigenvalues of L3 are denoted by = Ao < Ai < A2. 

2.2. Therefore Ai + A2 = 2(oi + 02 + 03). We claim that 01+02+03 > \/3 and 
the equality holds if and only if 6*1 = 6*2 = ^3 = f . 

Consider / := oi(0i) + 02(^2) + 03(^3) as a function defined on the domain 



^3 ■■= {{61,62, 63) \ei + d2 + e3 = TT,6^>o,i^ 1,2, 3}. 

To find the absolute minimum of /, we apply the method of Lagrange multiplier. 
Let 

F = 01 (0i) + 02(^2) + 03(^3) + yiOi +62 + 63- tt). 



Since 

dai{di) 



atli sm 6i 

we have 

dF 2^ 
0^^^-(l + a?) + y, 

0^-^-il + al) + y, 

0=§-^-il + al) + y, 
OF 

0= — =61 +62+63 -TT. 

dy 

Therefore the function / has the unique critical point {61, 62, 63) = (f, f , f )• 
Next, we investigate the behavior of the function / when the variable {61,62, 63) 

approaches the boundary of the domain ^3. Let {6i{t),62{t),63{t)),t € [0,cx)), be 

a path in the domain 573. Let ai{t) = cot 6i{t) for i = 1, 2, 3. 
Without loss of generality, we assume 

lim {6i{t), 62(1), 63(1)) = {0,S2, S3) 
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where S2 > 0, S3 > and S2 + S3 = tt. Then hmt^oo ai(t) — 00 and a2{t) + a3(t) > 
for t E [0, 00). Hence 

Um (ai(<) + a2{t) + a3{t)) = 00. 

If 9i + 6j < IT, then cot 9i > cot(7r — 9j) = — cot 9j. Therefore + aj > 0. Hence 
2/ = (ai + 02) + (a2 + 03) + (03 + ai) > 0. Thus / has the absolute minimum. But 
the absolute minimum can not be achieved at a point in the boundary of It 
must be achieved at the unique critical point ^, -1). 

This shows that 01+02+03 > -s/S and the equality holds if and only if di = 
02 = = ^. 



2.3. Since 

A2 = fli + 02 + 03 + V (oi + 02 + as)^ - 3, 

we have A2 > -\/3 and the equality holds if and only if 6*1 = 6*2 = ^3 = f . 

Since A1A2 = 3, we have Ai < a/S and the equality holds if and only if 9i ^ 62 = 

£^3 - 3- 



3. QUADRILATERALS 

3.1. The vertices of a cyclic quadrilateral decompose its circumcircle into four arcs. 
We assume the radius of the circumcircle is 1 and the lengths of the four arcs are 
26*1, 26*2, 26I3, 26*4. Let := cot 6*^ for i = 1, 4. The condition 61 + 62 + 6*3 + 6*4 = tt 
implies 



(3) 



040203 + 010204 + 0103O4 + 02O304 = 04+02 + 03 + 04. 



There are two ways to decompose a cyclic quadrilateral in to a union of two 
triangles. The two ways produce the same discrete Laplace operator: 



/ 04 + O4 — Ol 

— O4 Ol + 02 

-02 
—04 



The characteristic polynomial of L4 is 





-02 

02 + 03 

-03 



-04 ^ 



-03 
03 + 04 / 



P4{x) = — 2(oi + O2 + O3 + 04)2;'' 

+ (3(oi02 + 02O3 + 03O4 + 04O1) + 4(0403 + 0204))a;^ 



— 4(010203 + 01O2O4 + 01O3O4 + 020304)0; 



2(oi + 02 + 03 + 04).!;' 



+ (3(oi02 + 02O3 + 03O4 + 0401) + 4(0403 + 0204))a;^ 

— 4(oi + 02 + 03 + 04)0;, 



by the equation ([3]). 
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3.2. By the similar argument in the case of triangles, we can show that ai+a2+a3+ 
tti has the unique critical point at ^, f , ^). And we have 2(ai + 02 + 03 + 04) = 
(ai + 02) + (ai + as) + (03 + a^) + (04 + ai) > 0. 

Next, we investigate the behavior of the function oi + 02 + 03 + 04 when the 
variable {61,62, 03, 64) approaches the boundary of the domain 

f^4 = {(^1,^2,^3,^4) I ^1+^2 + ^3 + ^4 = 77,^^ >0,i = l,2, 3,4}. 

Let (0i(t),92(t),9a{t),94{t)), t e [0,oo), be a path in the domain O4. Let ai{t) = 
cot 6*1 (i) for i = 1,2,3,4. 

Without loss of generality, we assume 

lim {ei{t),e2{t),ei{t),ei{t)) = (o, S2,S3,S4) 

t— >-oo 

where Sj > for i = 2, 3, 4 and S2 + S3 + S4 = tt. And we can assume that S2 < f 
and S3 < ^. Then \imt^ooa.i{t) = 00, a2{t) > when t is sufficiently large and 
a3(t) + a4,(t) > for any t € [0, 00). Hence 

lim (aiCt) + 02 (t) + asft) + aM)) =00 

t— >cx> 

Therefore ai + 02 + 03 + 04 achieves its absolute minimum at the unique critical 
point (f , f , f , f )• Hence ai + 02 + 03 + 04 > 4 and the equality holds if and only 

if 0^ = 0^ = 03=04 = ^. 

Therefore Ai + A2 + A3 > 8, A1A2A3 > 16 and the equality holds if and only if 
0-f^ = 02 = 0^ = O4 = J. 

3.3. To verify the statement about A1A2 + A2A3 + A3A1, by the formula of the 
characteristic polynomial ^-4(0;), it is enough to show 

3(aia2 + 0203 + 0304 + 0401) + 4(aia3 + 0204) > 20 

and the equality holds if and only if = ^2 = ^3 = ^4 = f • 

In fact, consider g := 3(aia2 + 0203 + 0304 + 0401) + 4(aia3 + 02^4) as a function 
defined on the domain 0,4. 

To find the absolute minimum of g, we apply the method of Lagrange multiplier. 
Let 

G = 3(aia2 + 0203 + 0-304 + 0401) + 4(aia3 + 0204) + y{6i + 62 + O3 + O4 - tt). 
Then 

dG 

= ^ = -(3a2 + 3a4 + 4a3)(l + a?) + y, 
dG 

= — = -(3ai + 3a3 + 4a4)(l + al) + y, 
0V2 

dG 

= ^ = -(3a2 + 3a4 + 4oi)(l + al) + y, 
003 

dG 

= ^ = -(3a3 + 3ai + 4a2)(l + af) + y, 

UU4 

dG 

= — = ^1 + ^2 + ^3 + ^4 - TT. 

dy 

The first and the third equation above imply that 

(3a2 + 3o4 + 403) (1 + a\) = (3a2 + 3a4 + 4ai)(l + al) 
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which is equivalent to 

(ai — a3)(3aia2 + 3aia4 + 80203 + 80304 + 40103 — 4) = 0. 

We claim that the second factor is positive, i.e., 80102 + 801O4 + 80203 + 80304 + 
40103 > 4. 

In fact, since ^i + ^2 + ^3 < tt, then cot(^i + ^2) > cot(7r — ^3). Then 

0102 — 1 



oi + 02 



> -as 



which is equivalent to 

(4) 0102 + 0203 + 0301 > 1 

since oi + 02 > 0. 

By the similar reason, 

0104 + 0403 + 0301 > 1. 
At least one of 02 and 04 is positive. If 02 > 0, then 
80102 + 801O4 + 80203 + 803O4 + 4ai03 

= 8(0104 + 0403 + 0301) + (0102 + 0203 + 03O1) + 2(ai + 03)02 

> 8 + 1 + 0. 
If 04 > 0, then 

80102 + 801O4 + 80203 + 803O4 + 4ai03 

= (0104 + 0403 + 0301) + 8(0102 + 0203 + 0301) + 2(ai + 03)04 

> 1 + 3 + 0. 

Thus the only possibility is oi = 03 which implies 9i = 6*3. By the similar 
argument, = ^ and = ^ imply 62 = O4,. Since ^1+^2 + ^3 + ^4 = ti", we have 
Oi + O2 = which implies 01O2 = 1. 

Now = 1^ and = II imply 

(802 + 804 + 4a3)(l + Oi) = (801 + 803 + 4a4)(l + 02). 
Since Oi = 03 and 02 = 04, we have 

(602 + 4oi)(l + al) = (601 + 4o2)(l + oi). 
Since 01O2 = 1, we have 

(oi - a2)(a? +al + O1O2 + 1) = 0. 
Since the second factor satisfies 

aj + al + oioa + 1 = ^(a? + ^2) + ^(^i + (12)^ + 1 > 0, 

the only possibility is oi = 02. 

Therefore the function g = 3(ai02 + O2O3 + 03O4 + 04O1) + 4(ai03 + O2O4) has 
the unique critical point (f , f , f ■ f )• 

Next, we claim that 5 > 0. Since at least three of 01,02,03,04 are positive, 
without loss of generality, we may assume that Oi > 0, 02 > 0, 03 > 0. Let's write 

g = 2(0102 + 0204 + 0401) + 2(0203 + 0304 + 0402) + (02 + 04)01 + (01+04)03 + 40103. 

Then each term of sum above is positive. 
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At last, we investigate the behavior of g when the variable approaches the bound- 
ary of the domain n^. Let {6i{t),02{t),9^{t)^9i{t))^ t E [0, oo), be a path in the 
domain ^4. Let ai{t) = cot6'i(t) for i = 1,2,3,4. 

Without loss of generality, we have 

lim {9i{t),e2{t),e3{t),04{t)) = (0,S2,S3,S4) 

t^oo 

where Si > for i — 2/6,4 and S2 + S3 + S4 = tt. And we can assume that S2 < ^ 
and S3 < -I- 
Let's write 

g = 2{ai{t)a2{t) + a2{t)ai{t) + ai{t)ai{t)) 

+ 2{a2{t)a3{t) + a3{t)ai{t) + a4{t)a2{t)) 

+ {a2{t) + a4{t))ai{t) + {ai{t) + 04^)03^ + 4ai{t)a3{t). 

By the inequality (|3]), 

ai{t)a2{t) + a2{t)a4{t) + a4{t)ai{t) > 1 

and 

02(^)03(0 + a3{t)a4{t) + a4,{t)a2{t) > 1 
for any t G [0, 00). Since a2{t)+a4{t) > when t is sufficiently large, limt_^oo(a2(i) + 
a4{t))ai{t) = 00. And {ai{t) + a4{t))a3{t) > 0, 4ai(t)a3(t) > when t is sufficiently 
large. Hence g approaches 00. 

Therefore g has a lower bound and can not achieve its absolute minimum at a 
boundary point. It much achieve its absolute minimum at the unique critical point 

/ 7r 7r 7r tt \ 
"> 4 ' 4 ' 4 ' 4 /■ 

3.4. We verify the statement about Ai in this subsection. First, we verify that 
Ai < 2 as follows. Let 

P4(x) 

Q(x) := — ^ =x^ - 2{ai + 03 + + a4)x'^ 

X 

+ (3(aia2 + 0203 + a3a4 + a4ai) + 4(aia3 + a2a4))a; 

— 4(ai + 02 + 03 + 04). 

We have Q{0) = -4(ai + 02 + 03 + 04) < -16. 

If Q{2) > 0, then the first root of Q{x) is less that 2, i.e., Ai < 2. 

If (5(2) < 0, we claim that Q'(0) > and Q'{2) < 0. Once the two statements 
are established, Ai < A2 < 2. 

In fact Q'{0) — 3(aia2 + a2a3 + a3a4 + 0401) + 4(aia3 + a2a4) > 20. 

To verify Q'{2) < 2, we need to use the assumption Q{2) < 2. In fact Q{2) < 2 
implies 

3(aia2 + a2a3 + a3a4 + a4ai) + 4(aia3 + 0204) < 6(ai + 02 + a3 + 04) — 4. 
Now 
Q'(2) 

= 12 — 8{ai + 02 + a3 + a4) + 3(aia2 + a2a3 + a3a4 + a4ai) + 4(oia3 + 0204) 
< 12 — 8(ai + 02 + 03 + 04) + 6(ai + 02 + 03 + 04) — 4 
= 8 — 2(ai + a2 + 03 + 04) 
<0, 
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since ai + a2 + 03 + 04 > 4. 

Second, wc verify that Ai = 2 if and only ii 61 = 62 = 63 = 04 = j. Since Ai = 2 
is the first root of Q{x), we have Q'{2) > 0. On the other hand, it is shown that 
(5(2) < imphes Q'{2) < 0. Hence the only possibility is Q'{2) = 0. This requires 
that ai + a2 + 03 + 04 = 4. Therefore we must have Oi = 62 = O3 = 64 = 

4. GENERAL CYCLIC POLYGONS 

4.1. Assume n > 5. The vertices of a cyclic n-gon decompose its circumcircle into 
n arcs. We assume the radius of the circumcircle is 1 and the lengths of the n arcs 

are 2^1, 26*2, 26*,,,. 

The discrete Laplace operator of a cyclic n-gon is independent of the choice of a 
triangulation. It is 



ai + a„ 


-ai 








-a 


-ai 


ai + 02 


-02 











-02 


012 + as 


-as 











-03 


03 + 04 














—04 






V 







The eigenvalues are = Aq < Ai < 




• • < A„_i. 



-1 + On / 



4.2. Wc have Y^"^i = 2^^^^ai. By the similar argument in the case of tri- 
angles and cyclic quadrilaterals, we can show that X)r=i has the unique critical 
point (01,..., 0„) = (^,...,^). 

Since there is at most one non-positive number in ai,...,a„, without loss of 
generality, wc mav assiunc 04 > 0, ...,a„_i > 0. Since a„_i -|- a„ > 0, wc have 

Er=i««>o. 

We investigate the behavior of Y^"^i ai when the variable approaches the bound- 
ary of the domain 



{{9i,...,0„) I ei + ... + en = 7T,ei>o,t = i,...,n}. 



Let {9i{t),92{t),9z{t),0{t)), t G [0,00), be a path in the domain ^4. Let ai{t) = 
cot 6*1 (t) for i = 1,2,3,4. 

Without loss of generality, we have 

lim {0i{t),...,9n{t)) = {Q , S2) "M ^n) 

t—>oo 

where Sj > for i = 2,...,n and S2 + ... + s„ = tt. And we can assume that 
S2 < f , Sn-i < f • Since aj(i) > for i = 2, n — 1 and a„_i -|- a„ > when t 
is sufficiently large, limt_>.oo ai = 00 implies that limt^oo J27=i ^« ~ 
Thus Yl"--, tti achieved the absolute minimum at (J, J). 

4.3. In this subsection we verify the statement about 11^=/ 

The weighted matrix-tree Theorem. Let M be an n by n matrix. If the sum 

of the entries of each row or each columm, of M vanishes, all principle n — 1 by n—1 
submatrices of M have the same determinant, and this value is equal to ^ times 
the product of all nonzero eigenvalues of M. 
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For the reference of the weighted matrix-tree Theorem, for example, see |LWj . 
page 450, Problem 34A or |DKMj . Theorem 1.2. 

In our case, according the weighted matrix-tree Theorem, to calculate 0"=/ 
of the matrix L„, it is enough to calculate a particular principle n — 1 by n — 1 
matrix. 

Lemma 4. Let Nn be the submatrix obtained by deleting the first row and first 
column of the matrix Ln ■ Then 

n 

det A^„ = ai...ai...an, 

i=l 

where Oi means that Oi is missing. 

Proof. We prove the statement by the mathematical induction. It holds for n = 4 
as we see in the formula of the characteristic polynomial Pi{x). We assume it holds 
for n < m — 1. By the property of tridiagonal matrices, we have 

det N„r = (a„i_i + a„i) dot iY„i_i - al^_i det N„i-2- 

Then by the assumption of the induction, 

m— 1 m— 2 

detNm = (flm-i + flm) ^ fli . . .a^ . . .Om- 1 - a'^_^ ^ ai...ai...am-2 

i=l i=l 
m—1 ra—2 

= o-m-i ai...ai...am-i — am-i ai...ai...am-20,m-i 

i=l i=l 

m—1 

+ a?ri ai...ai...am-i 

1=1 

m—1 



ii...am-i + a 



ai...ai...a„ 

i=l 



□ 

In the following, we prove that X]"=i ai...ai...a„ achieves its absolute minimum 
when di = ... = 9n = ^. It is enough to show that 

a. X]"=i ai...ai...an has the unique critical point ( — ,...,—); 

Yl7=i ai...ai...an > 0; 
c. ^"^-^ ai...ai...a„ approaches oo as the variable approaches the boundary of 
the domain 

When n = 4, since 010203 + aia2a4 + aio^Oi + 020^04 — ai + 02 + 03 + 04, the 
three statements above are already shown to be true in section 3. We assume that 
the three statements above hold when n < m — 1. 

Let's check the three statements when n — m. Consider the function 

m 

H = ^oi. ..ai...am - v{Oi + ...Om - ti"). 

i=l 
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Then = ||- and = imply that 

m m 

{as. ..am + ^ as... Si... am) + a?) = {as. ..am + X] «3"-aj...am)(l + al). 

i=S i=3 

Since ai + 02 > 0, it is equivalent to 

I \/ \l CHCI2 — i \ ^ ^ \ 

(ai - ai){a\ -\- a2){as...am H ; > as...ai...am) = 0. 

ai + 0,1 

1=3 

The third factor is 

m 

cot ^3... cot dm + COt(0i + Q-2) ^^COt 03...COt ^j... COt 0^ 

i=3 

which is written as ^^^^^'a\..Mi..Mm-\i where a\ = cot(0i + 92),ai = cot^^+i for 
1 = 2, 1. This expression corresponds to a cyclic (to— l)-gon. By assumption 

of the induction, ^^'^-^ ai-.-aj-.-am-i > 0. 

Hence the only possibility is ai = ai. By similar argument, we show that ai = aj 
for any Hence the function ^ a\...ai...am has the unique critical point such 
that 0i = ^ for any i = 1, ...,m. 

Next, we claim that Y^^i ai---ai...am > 0. Without loss of generality, we assume 
that ai > 0, a2 > 0, am-i > 0. Now 

m 

^ai...ai...am 



i=l 

m-2 



m—Z 

= aia2...am-2{(lm-l + dm) + ^ ai...ai...am-2{am-l(lm) 

i=l 

m—2 m—2 

= aia2...am-2{am-i + am) + ^ ai...ai...am-2{am-iam - 1) + X] c^i-^i-^n 

i=l i=l 
m—2 ^ m—2 

= (a^-i + am){aia2...am-2 + ai...ai...am-2 ) + ^ ai...ai...a, 

1=1 am-i + am .^^ 

Let am-\ = "a'l^+arl = cot(6'm_i + Then 

m-2 . 
I \ ^ ^ <^m— l<^m ~ J- 

aia2...aTO-2 + ai...a8...aTO_2 — 

^^m— 1 ~r ^m 

m-2 

= aia2...am-2 + X^ ai...ai...aTO_2am-i- 
1=1 



Consider an cyclic (to — l)-gon with angles 0i, Om-2, Om-i + (^m- By the assump- 
tion of induction, 

m-2 

aia2...am-2 + X] ai— ai-am-2am-i > 0. 

i=l 

Therefore Yl^i a\...ai...am > 0. 
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At last, we investigate the behavior of the function X^I^i o,i---ai...am when the 
variable approaches the boundary of the domain 

= {{91, e,n) I 01 + ... +0m^TT, 9, > 0, I = 1, m}. 

Let (6*1 (t), f?,n(t)), t e [0,oo), be a path in the domain Q,m- Let ai{t) = cot6'i(t) 
for t = 1, TO. 

Without loss of generality, we assume 

lim {ei{t),...,er,i{t)) = (0, S2,...,Sm), 

where S2 > 0, Sm > and S2 + ... + Sm = tt. And we can assume furthermore 
that S2 < ^,...,Sm_i < ^. Thus ai(i) > 0, am-i(i) > when t is sufficiently 
large. To simplify the notation, we denote ai(t) by Ui in the follows. Now 

rn 

ai...ai...am 

4=1 

m — 2 m — 2 

= (a,„_i + a„i){aia2...am^2 + ^ ai...ai...am_2-^^^^ j ) + ^ ai...ai...am-2 

1=1 2=1 



m — 2 m — 2 



(a,„_i + a„i){aia2...am-2 + ^ ai...ari...am_2am_i) + ^ ai...a;;...a 



^ ^ m-2, 

i=l 

where a„i_i = = cot(6',„_i + e^). 

By the assumption of induction, 

m-2 

a\_a2--am-2 + ^ ai...ai...am_2am-i > 

for any t g [0, cxj). Since Cm-i + Qm > for any t G [0, oo) and X]"=i^ ai...ai...am-2 
approaches co, we see that X^IHi 0'\---o^i---o-m approaches oo. 
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